^ 



The A^-soliton solution of a two-component 
modified nonlinear Schrodinger equation 

Yoshimasa Matsunqf 
Division of Applied Mathematical Science, 
Graduate School of Science and Engineering 
P^ ■ Yamaguchi University, Ube, Yamaguchi 755-8611, Japan 

xn : ABSTRACT 

The iV-soliton solution is presented for a two-component modified nonlinear Schrodinger 
^ \ equation which describes the propagation of short pulses in birefringent optical fibers. The 

solution is found to be expressed in terms of determinants. The proof of the solution is 
carried out by means of an elementary theory of determinants. The generalization of 

the 2-component system to the multi-component system is discussed as well as a (2+1)- 

> ■ 

\0 ' dimensional nonlocal equation arising from its continuum limit. 

!>■ 

p 

O ; PACS: 05.45.Yv; 42.81. Dp; 02.30.Jr 

Keywords: Modified nonlinear Schrodinger equation; A^-soliton solution; Two-component 
system 



"^E-mail address: matsuno@yainaguchi-u.ac.jp 



1. Introduction 

We consider tlie following two- component system of nonlinear evolution equations 
which is a hybrid of the coupled nonlinear Schrodinger (NLS) equation and coupled deriva- 
tive NLS equation 

(2 \ r / ^ \ " 

= 0, (j = l,2), (1.1) 





1 9i,t + (lo^xx + /" >^ l^fc I ^i + 17 



where q^ = qj{x, t) {j = 1, 2) are the slowly varying complex envelopes for polarizations, /i 
and 7 are real constants and subscripts x and t appended to Qj denote partial differentia- 
tions. This system of equations has been derived as a model for describing the propagation 
of short pulses in birefringent optical fibers [1] . Its integrability has been established by 
constructing the Lax pair and an infinite number of conservation laws [2]. Hence, the 
initial value problem can be formulated in principle by means of the inverse scattering 
method. A few exact solutions are now available for the system (1.1). Actually, one- and 
two-soliton solutions have been presented by employing Hirota's bilinear transformation 
method [3, 4]. The A^-soliton solution of the system (1.1) with /i = was also obtained 
quite recently by means of the Darboux transformation where N is an arbitrary positive 
integer [5]. 

The purpose of this Letter is to construct the bright A^-soliton solution of the system 
(1.1) within the framework of the bilinear formalism. We present a compact determinantal 
expression of the A^-soliton solution which would be suitable for investigating the inter- 
action process of solitons. The proof of the solution is performed by using an elementary 
theory of determinants. In concluding remarks, the generalization of the system (1.1) to 
the n-component system {n > 3) is briefly discussed as well as a (2+l)-dimensional non- 
local modified NLS equation arising from its continuum limit. Furthermore, some issues 
associated with the system (1.1) will be addressed. 

2. Bilinear form and A^-soIiton solution 

2.1. Bilinear form 



We first apply the gauge transformations 



Qj = Uj exp 



^ -^-^ fc=i 



, (J = l,2), 



(2.1) 



to tlie system (1.1) and transform it to the system of equations for -Uj 



lUj^t + Mj>x + At I '^ l^fcl^ ) U'J + 17 I y^^fc^fc.a: ) ^j = 0, (j = l,2), (2.2) 

where the asterisk denotes complex conjugate. Note that we have imposed the boundary 
conditions q^ — )■ 0,-Uj — )■ as |a;| — )■ oo to give the bright soliton solutions. By means of 
the dependent variable transformations 



w, = y, (J = l,2), 



(2.3) 



the system (2.2) can be transformed to the following system of bilinear equations for / 
and Qj [4] 

(iA + B^ft ■ / = 0. (J = l,2), 



Dj-r^'^Y. 



9k\ 



k=l 



k=l fc=l 

Here, the bilinear operators D^ and Dt are defined by 



(2.4) 
(2.5) 

(2.6) 



where m and n are nonnegative integers. It follows from (2.3) and (2.5) that 



(2.7) 



i'/ V^ I ,2 ^1 /* 

— > Ml. = -— m — , 



17 
2 



(2.^ 



fc=i 



which, substituted into (2.1), yields the solution of the system (1.1) in the form Qj 
g,r/P{j = 1,2). 



2.2. N-soliton solution 



We now state the main result in this Letter: 

The N-soliton solution of the system of bilinear equations (2.4), (2.5) and (2.6) is 
given by the determinants f and Qj (j = 1,2) where 



f 



A I 
-I B 



9j 



A 







/ 
B 







(J = l,2). 



(2.9) 



Here, A, B and I are N x N matrices and z, a^ and are N -component row vectors 
defined below and the symbol T denotes the transpose: 



A 



B 



[0'jk)l<j,k<N, O-jk 



1 ZjZj^ 



* ' ^3 



2pj+pl 
{fi + i'jpk)cjk 



exp{pjX + ip]t), 



Pi+Pk 



(^jk 



E 



{bjk)l<j,k<N, bjk 

{Sjk)i<j,k<N, '■ N X N unit matrix, 





'^sj'^ski 



(2.10a) 

(2.106) 

(2.10c) 
(2.10rf) 



z = {zi,Z2,...,zn), a.j = {aji,aj2,...,ajN), = (0, 0, ..., 0). 

The above A^-sohton solution involves 3A^ complex parameters pj and aij,a2j {j = 

1.2, ...,N). The former parameters determine the amplitude and velocity of the solitons 
whereas the latter ones determine the polarizations and the envelope phases of the solitons. 
In the special case of A^ = 1,2, (2.9) and (2.10) reproduce the 1- and 2-soliton solutions 
presented in [4]. Note, however that the compact determinantal expression of the solution 
has been obtained here for the first time. It is also remarked that in the single-component 
case (i.e., the system (1.1) with g2 = 0), our solutions reduce to the soliton solutions given 
in [6]. 

2.3. Remark 

If we introduce the transformations 



/ = /, 9j = exp 




93, (J = l,2) 



2/i ~ 

•Jy Jy I (/ • (j (/ • 
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(2.11a) 
(2.116) 



then the bihnear equations (2.4), (2.5) and (2.6) recast to 

{iDi + Dlrg,-f = 0, (j = l,2), (2.12) 

2 

DJ-f* = jY.\9k\', (2.13) 

fc=i 



2 



Dlf-r = ^J2^^9k-9l (2.14) 

fc=i 

respectively. Thus, the form of Eqs. (2.4) and (2.5) is invariant whereas (2.6) becomes a 
simphfied equation with n = 0. Consequently, the proof of the A^-soliton solution may be 
performed for the corresponding solution with /i = 0. This reduces the total amount of 
calculations considerably, in particular for Eq. (2.6). Hence, in the analysis developed in 
Sec. 3, we put /x = without loss of generality. 

3. Notation and formulas for determinants 

To perform the proof of the A^-soliton solution concisely, we introduce some notations 
and write the basic formulas for determinants. 
3.1. Notation 

We define the following matrices associated with the iV-soliton solution (2.9): 

D - (_^, ^) , (3.1) 

/A I 0^\ 
D(a;b)= -/ B b^ , (3.2) 



A 


I\ 




-I 


b) 


f 


A 


I 


(f' 


-I 


B 


b^ 





a 


^ 


A 


/ 


z^^ 


-I 


B 


0^ 





a 


0. 


A 


/ 


z^ 


-I 


B 


0^ 


z* 









L'(a;z)= -/ B 0^ , (3.3) 



D(z*;z)= -/ B (f\. (3.4) 

\z* 0/ 

Note the position of the vectors a, b, z and z* in the above expressions where a and b 

represent a* and a^ (j = 1,2), respectively. The matrices which include more than two 

vectors will be introduced as well. 



3.2. Formulas for determinants 

In the process of the proof, we use some basic formulas in addition to the fundamental 
properties of determinants, which we shall now write down for later convenience. Let 
^ = {0'jk)i<j,k<M be an M X M matrix with M being an arbitrary positive integer and 
Ajk be the cofactor of the element ajk- Then, we have [7] 



(3.5) 



(3.6) 

|A(ai,a2;bi,b2)||A| = |A(ai; bi)||A(a2; ba)! - |A(ai; b2)||A(a2; bi)|. (3.7) 

The formula (3.5) is the differentiation rule of the determinant and (3.6) is the expansion 
formula for a bordered determinant with respect to the last row and column. The formula 
(3.7) is Jacobi's identity and it will play a central role in the proof. 
The following two formulas are also employed frequently: 

|A(ai;bi)| ■■• |A(ai;b„; 



d 

dx 


j,k=l 


A a^ 
b z 


M 

= A\z- y ' Ajkttj 
j,k=l 



|A(ai,...,a„;bi,...,b„)||A| 



n-l 



|A(a„;bi) 



|^(a„;b„) 



(n>2) 



(S.t 



|A + ^bfa,| = |A|-5^|A(a,;b,)| + |A(ai,a2;bi,b2)|. (3.9) 

s=l s=l 

The formula (3.8) can be proved by a mathematical induction. The case n = 2 reduces 
to Jacobi's identity (3.7). The notation bj'ag on the left-hand side of (3.9) represents an 
M X M matrix whose (j, k) element is given by hgjask- The formula (3.9) is derived simply 
using the property of the bordered determinant. 

4. Proof of the A^-soliton solution 

In this section, we perform the proof of the A^-soliton solution by means of an ele- 
mentary theory of determinants. To this end, we need the differentiation rules of the 
determinants / and Qj with respect to x and t as well as their complex conjugate expres- 
sions. We first write down them and then proceed to the proof. 



4-1- Differentiation rules of f and Qj 

In terms of the notation introduced in Sec. 3.1, we have 

f=\D\, g, = -\D{8i*;z)\. (4.1) 

The differentiation rules for these determinants are derived simply with use of the formulas 
(3.5) and (3.6). They read 

/* = -^{|^(z*;z.)|-|Z^(z:;z)|}, (4.2) 

/, = -i|Z}(z*;z)|, (4.3) 

/.x = -^{|^(z*;z.)| + |Z^(z:;z)|}, (4.4) 

g,,t = -\D{s.*;zt)\ + '-\Dia.*,z*;z,z^)\, (4.5) 

^j,, = -|D(a*;z,)|, (4.6) 

9j,xx = -|I^(a*;z^^)| + -|D(a*,z*;z^,z)|. (4.7) 

For an instructive purpose, we give a proof of (4.3). The other formulas can be proved in 
the same way. First, let D = {djk)i<j,k<2N and Djk be the cofactor of djk- We then apply 
the formula (3.5) to the determinant / given in (2.9) to obtain /^ = | Xli a:=i -^ife-^i-^fc- 
This can be put into the form (4.3) in view of the formula (3.6) with z = and the 
definition (3.4). The key in the proof is that the factor {pj + pl)^^ in the element ajk of 
the matrix A has been cancelled after differentiation with respect to x. 

We compute the complex conjugate expressions of these differentiations. To this end, 
we note from (2.10) that A* = A^ and B* = B^ — i^yC'^ where C is an A^ x A^ matrix 
with elements Cjk defined by (2.106). These relations lead to the expression of /* 



f* 



A I 

-I B-i-fC 



{U 



Using the formula (3.9), we can expand /* in powers of 7 



2 



r = |D|+i7}^|D(a*;a,)|-7'l^(aI,a;;ai,a2)|. (4.9) 



The similar procedure applied to f* and g* gives 

-1-2 2 

/: = --|D(z*;z)|-^^|D(z*,a*;z,a,.)l + ^l^(z*,at,a;;z,ai,a2)|, (4.10) 



i=i 



^* = \D{z*;aj)\ +i7^|i:'(z*,a^;aj,afc)|. 



(4.11) 



fc=i 



The above formulas will be used to prove Eqs. (2.4) and (2.5). For the proof of 
Eq. (2.6), on the other hand, we need another differentiation rules, which we shall now 
describe. First, we rewrite / in the form 



/ 



A I 
-I B 



where A and B are N x N matrices defined by 

A = {djk)i<j,k<N, o,jk 

B = {hjk)i<j,k<N-, bjk = 



1 1 



hCjkPk * 

Pi+Pk ^ 



Using the formulas (3.5) and (3.6), we obtain an alternative expression of fx 



(4.12a) 



(4.126) 
(4.12c) 



/, = -i7 5^|/^(b*;a,)|, 
i=i 

where bj is an A^-component row vector given by 

hj = {ajipl, aj2P*2, ..., ajNP*N)- 
Differentiation of (4.13a) with respect to x gives 



(4.13a) 



(4.136) 



fxx = —'^\D{z*,h*;z,aj)\, 



Similarly, we have 



(4.14) 



9j,x 



\D(h*; z)| + i7 ^ |D(a*, b*; z, a^) 



k=l 



(4.15) 



The complex conjugate expressions of /^, f^^ and Qj^^ are derived by referring to the 
procedure used to obtain /*. They read 

2 2 

/: = i7j]|I?(a*;b,)|-7^ J]|D(a*,a-b,,afc)|, (4.16) 

j=l j,k=l 

2 2 2 

fxx = -Y^\Diz*,a*;z,hj)\ + YYl l^(z*' ^i' ^^5 z, bj, afe)|, (4.17) 

i=l j,k=l 

2 

9j,x = \D{z*;hj)\+i-f'^^\D{z*,cil;hj,ak)\ + |D(z*, a^; a^, bfc)| j 

fc=i 

2 

'||D(z*,a*,a;;b,-,ai,a2)|+ J]] |/^(z*, a*, a^; a,-, b^, a^l}. (4.18) 



fc,Z=l 

It follows from (4.3) and (4.13a) that 



|D(z*;z)| = 2i7 5^|I^(b*;a,)|. (4.19) 

i=i 

Equating (4.10) and (4.16) and using (4.19) to eliminate \D{z*;z)\, we find that 

2 2 

J]{|D(a*;b,)| + |/^(b*;a,)|} = --5^|/^(z*,a*;z,a,)| 

2 .. 

-i7| J2 |^(a*,a*;bj-,afe)| + -|D(z*,at,a;;z,ai,a2)|}. (4.20) 

j,k=i 

4.2. Proof of Eq. (2.4) 

Let Pi be 

Pi = (iA + Dl)g, ■ f. (4.21) 

Substituting (4.1)-(4.7) into (4.21), Pi becomes 

Pi = -|Z}(a*,z*;z,z,)||Z}| + \D{si*;z)\\D{z*;z,)\ - \D{si*;z,)\\D{z*;z)\ 

-{i|D(a*;Zi)| + |/^(a*;z,,)|}. (4.22) 



If we use Jacobi's identity (3.7) and the fundamental formula a|D(a; bi)| +/3|D(a; b2)| = 
|D(a; ahi + /3h2)\ (a,/3 G C), Pi simplifies to 

Pi = -|D(a*;iZi + z,,)|. (4.23) 

Since iz^ + z^x = by (2.10a), the last column of the determinant consists only of zero 
elements. Hence, we see that Pi = 0, which completes the proof of Eq. (2.4). 

4.3. Proof of Eq. (2.5) 

Let P2 be 

2 

fc=i 
If we substitute (4.1), (4.3) and (4.9)-(4.11) into (4.24) and use Jacobi's identity (3.7), we 
can reduce P2 to 

P2 = y{ - |I)(z*,at,a;;z,ai,a2)||Z^| + |D(z*; z)||D(at, a^; ai, aa)] 

-|D(at;z)||D(z*,a;;ai,a2)| + |/^(a;;z)||Z}(z*,at;ai,a2)|}. (4.25) 

Assume that the matrix D is non-singular, i.e., |Z^| 7^ which ensures the regularity of 
Uj from (2.3). Multiplying \D\ on both sides of (4.25) and using Jacobi's identity (3.7) as 
well as the formula (3.8) with n = 3, we see that |P'|P2 = 0, implying that P2 = 0. This 
completes the proof of Eq. (2.5). 

4.4. Proof of Eq. (2.6) 

The proof of Eq. (2.6) is very complicated compared with that of Eq. (2.4) and (2.5) 
although it can be performed straightforwardly. Let P3 be 

2 

P, = Dlf.f*-^Y.^^9k-9l. (4.26) 

fc=i 

By using (4.1) and (4.13)-(4.18) and applying Jacobi's identity (3.7) and the formula (3.8) 

with n = 3, P3 can be written, after some lengthy calculations, in the form 

2 • 3 

Pz = jPsi + yP32 + ^Pss, (4.27a) 
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with 



P3i = |/^(z*;z)|^{|Z}(a*;b,)| + |D(b*;a,)|}, (4.276) 

2 

A2 = -$^{|/^(b*;a,)||/^(z*,a*;z,a,)| + |D(z*;z)||D(a*,a*;b„afe)|}, (4.27c) 
j,k=i 

2 

^33 = X^{ - l^(at,a;;ai,a2)||D(z*,b*;z,aj)| - |L'(z*, a^; ai, a2)||L'(a;, b*; z, a^-)! 



+ |L'(z*, al; ai, as)] |L'(ai, b*; z, a^) 
If we introduce (4.19) and (4.20) into P31, P3 becomes 

„3 



(4.27rf) 



• 3 -^ 

P3 = — ^||P'(b*;aj)||P'(z*,at,a;;z,ai,a2)| - |P'(a*, a^; ai, a2)||P'(z*, b*; z, a^ 



~\D{z*, a*; ai, as)! |Z^(a;, b*; z, a^OI + l^(z*, a^; ai, as)! |Z^(a*, b*; z, a^-)! }• (4-28) 

We can confirm that !-D!^P3 = by virtue of Jacobi's identity (3.7) and the formula (3.8) 
with n = 3. Thus, we complete the proof of Eq. (2.6). 

5. Concluding remarks 

5.1. Generalization to the n-component system 

The generalization of the system (1.1) to the n-component system {n > 3) will be 
obvious. It takes the form 



i <lj,t + qj,xx + f^ \^\(lk\^ ) Qj + 17 



vfe=i 






Qk\ Qj 



0, (j = l,2,...,n). (5.1) 



By means of the gauge transformation similar to (2.1), one can transform (5.1) to the 
system (2.2) in which the summation with respect to k extends from 1 to n. The system 
of bilinear equations corresponding to (2.4), (2.5) and (2.6) can be derived easily. They 
read 



{iD, + Dl)gyf = 0, (j = l,2,.... 



n 



Dj-r = 'iYl 



\9k\ 



(5.2) 
(5.3) 



fc=i 
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n . n 



Dlf ■ r = /i^ l^^l' + 'iJ2^-9k ■ 91 (5.4) 

k=l k=l 

The iV-soliton solution of the above system of bihnear equations would be given by (2.9) 
with the element Cjk in (2.10b) being replaced by X]r=i '^^i'^sfc- The key formulas in 
performing the proof of the A^-soliton solution are the formula (3.8) and the generalization 
of the formula (3.9) 

n n' 

l^ + $^bfa,| = |A| + ^(-l)™ Y, |A(a,,,a,„...,a,„;b,,,b,,,...,b,^)|, 

S=l m=l l<Si<S2<...<Sm<™ 

(5.5) 
where n' = min(n, M). 

5.2. Nonlocal modified NLS equation 

Another interesting issue to be worth studying is a nonlocal coupled NLS and deriva- 
tive NLS equation arising from a continuum limit of the system (5.1), which is written 
as 

iqt + Qxx + f^i \q\'^dz] q + i'y i \q\^dzq] =0, q = q{x,t,z). (5.6) 

When 7 = 0, this equation coincides with a (2+l)-dimensional nonlocal NLS equation 
proposed by Zakharov [8]. Its A^-soliton solution has been obtained in [9]. The bilinear 
equations corresponding to (5.2), (5.3) and (5.4) stem from their continuum versions where 
gj is replaced hj g = g{x, t, z) and the sum with respect to k chages to the integral with 
respect to z. We can expect that the A^-soliton solution of Eq. (5.6) takes the same form as 
(2.9) with Cjk in (2.10b) being replaced by j'^^aj{z)a*f,{z)dz where aj{z) (j = 1,2, ...,n) 
are continuous functions of z. However, its proof still remains open. The detailed analysis 
of the problems mentioned above will be done in future works. 

5.3. Miscellaneous remarks 

1. The results in this Letter have wide applications in the description of the propagation 
characteristics of ultrashort pulses in an optical fiber. One of which has been presented 
in [4] for investigating the dynamics of two optical solitons. The interaction process of an 
arbitrary number of solitons would be worth studying. As for various applications, one 
can refer to the references [1], [3] and [4] as well as a number of papers cited therein. 
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2. The extended or higher-order NLS equations have been studied previously by means 
of the bihnear method. Unhke the system (1.1), they require the third-order dispersion. 
On can consuh the references [10] and [11], for example as for the construction of the 
A^-soliton solutions . 

3. An advantage of using the determinantal expressions of the solution is that the proof of 
bilinear equations reduce to the well-known formulas for determinants, as demonstrated 
in this Letter. On the other hand, it is possible to represent the solution in the form of 
a finite sum of exponential functions. See, for example, an A^-soliton solution of a single- 
component modified NLS equation in which the proof of the A^-soliton solution has been 
performed by an ingenious mathematical induction [6]. Its proof is, however, not easy to 
understand when compared with an elementary proof relying on Jacobi's identity. 

4. The repeated application of the Darboux transformation to the system (1.1) with 
/i = yields the bright A^-soliton solution of the form qi = {H2/ Hi)^, q2 = {H^/Hi)^ 
where Hj{j = 1,2,3) are 3A^ x 3A^ determinants [5]. In terms of the bilinear operator, 
they can be rewritten as gi = DxH2 ■ Hi/Hf,q2 = D^H^ ■ Hi/Hf. Comparing these 
solutions with the A^-soliton solutions (2.9) with /x = 0, one can infer the relations / = 
cHi, f*gi = c^DxH2 ■ Hi and f*g2 = c^D^H^ ■ Hi with c being a complex constant. The 
proof of these relations seems to be a highly nontrivial problem to be pursued in a separate 
context. 
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